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ABSTRACT 

We examine the relationship between source radio luminosity and kinetic power in Active Galactic 
Nucleus (AGN) jets. We show that neglecting various loss processes can introduce a systematic bias 
in the jet powers inferred from radio luminosities for a sample of radio galaxies. This bias can be 
corrected for by considering source size as well as radio luminosity; effectively the source size acts as 
a proxy for source age. Based on a sample of FR-II radio sources with jet powers derived from the 
measured hotspot parameters, we empirically determine a new expression for jet power that accounts 



for the source size, 
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(l + ^rfe) 



0.58±0.17 



where D is source size 



lO^eW ~ -^-^-0.8 llO^^WHz-i, 

and ii5i the 151 MHz radio luminosity. By comparing a flux-limited and volume-limited sample, we 
show that any derived radio luminosity - jet power relation depends sensitively on sample properties, 
in particular the source size distribution and the size-luminosity correlation inherent in the sample. 
Such bias will affect the accuracy of the kinetic luminosity function derived from lobe radio luminosities 
and should be treated with caution. 
Subject headings: galaxies: active — galaxies: jets — quasars: general 



1. INTRODUCTION 

There has been much recent interest in accurate mea- 
surement of the kinetic power of Active Galactic Nu- 
cleus (AGN) jets. Such measurements facilitate studies 
of AGN physics (for example jet composition, and gen- 
eration and collimation mechanisms); as well as helping 
to quantify the effects of AGN feedback on galaxies and 
their building blocks in a cosmological framework. 

AGN jet powers in individual objects have been esti- 
mated using cavities in the hot X-ray gas (RafFerty et al. 
2006; Birzan et al. 2008; Cavagnolo et al. 2010); dynam- 
ical models of radio sources (Machalski et al. 2004; Sha- 
bala et al. 2008); and radio lobe expansion speeds (Daly 
et al. 2012). Recently, in a companion paper (Godfrey 
& Shabala 2012, hereafter Paper 1) we presented a new 
method based on measured properties of jet termination 
shocks known as hotspots. This method agrees well with 
other techniques. 

It is highly desirable to develop a very simple method 
of measuring AGN jet power, ideally one that relies on a 
minimal number of easily measurable observables. The 
AGN radio luminosity function (RLE) is the standard 
measurement used to study cosmological evolution of the 
AGN population. Various authors (e.g. Willott et al. 
1999; Heinz & Sunyaev 2003; Merloni & Heinz 2007; 
Cavagnolo et al. 2010) have attempted to convert radio 
luminosity, Lradio, to AGN mechanical power, Qjct. Such 
a conversion is not easy, since the observed radio lumi- 
nosity depends on the details of the electron acceleration 
and loss processes, as well as interactions between the 
radio source and its environment. However, the stakes 
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are high, since a correct iradio — Qjet relationship allows 
the whole radio luminosity function to be converted to 
a kinetic jet power function (e.g. Merloni & Heinz 2007; 
Martinez-Sansigre & Rawlings 2011; Kaiser & Best 2007; 
Cattaneo & Best 2009), and therefore AGN heating rate. 
Cattaneo & Best (2009) show that the uncertainty in the 
assumed Lradio — Qjet relation dominates the uncertainty 
in the kinetic jet power function. It is therefore of funda- 
mental importance to determine the -/jradio~Qjot relation, 
or relations. The ultimate goal is to self-consistently rec- 
oncile in the framework of galaxy formation models the 
observed AGN properties (i.e. the RLE) with galaxy and 
interstellar/intracluster medium properties (mass func- 
tions, colours). 

One of the most commonly used relations in converting 
radio luminosity to jet power is that presented by Willott 
et al. (1999). Based on synchrotron minimum energy 
calculations in combination with a self-similar model of 
radio galaxy evolution (Falle 1991; Kaiser & Alexander 

1997), these authors predicted Qjet oc i^adio- This pro- 
portionality was supported by an observed correlation 
between [O II] narrow-line luminosity, assumed to be a 
proxy for the continuum ionising luminosity and there- 
fore related to jet poweiQ, and 151 MHz lobe radio lu- 
minosity. Unlike many other correlations (e.g. Birzan 
et al. 2008; Daly et al. 2012) this result is not purely 
empirical: the correlation can be understood in terms of 



^ As we point out in Paper 1, the slope of the Lqii ~ ^151 corre- 
lation is highly uncertain and sample dependent (e.g. Hardcastle 
et al. 2009, Fcrnandes et al. 2011). Furthermore, how well the 
narrow-line luminosity traces jet power is highly debatable (com- 
pare e.g. Hardcastle et al. (2009) and Tadhunter et al. (1998) with 
Rawlings & Saunders (1991) and Shabala ct al. (2012)). 
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a dynamical radio source model under some quite rea- 
sonable assumptions. 

In this work, we use the hotspot method described in 
Paper 1 to re-examine the Qjot — ^radio relation of Willott 
et al. (1999). We show that the conversion of radio lu- 
minosity to jet power is sensitive to the size of the AGN 
radio lobes, largely due to inverse Compton losses. We 
note that the analysis of Willott et al. (1999) neglected 
this important radiative loss mechanism, which becomes 
increasingly relevant for old (> 100 Myrs) AGN, result- 
ing in a systematic shift in jet power estimates. As a 
result of these losses, the slope and normalisation of the 
Qjet — iradio relation are dependent on the distribution 
of source sizes within the sample being considered. 

Throughout the paper, we adopt the concordance cos- 
mology of nu = 0.27, r^A = 0.73, h = 0.71. 

2. JET POWER FROM SOURCE DYNAMICS 

2.1. Radio source dynamics 

The class of dynamical models of double-lobed 
Fanaroff-Riley Type II radio sources, pioneered by 
Scheuer (1974), and subsequently developed by Begel- 
man & CiofB (1989), Falle (1991) and Kaiser & Alexan- 
der (1997) among others, describe the evolution of 
an overprcssured cocoon of radio plasma. Jets of 
synchrotron-emitting plasma emanating from the black 
hole accretion disk smash into the ambient gas form- 
ing termination shocks; observationally these are man- 
ifested as bright regions of radio emission, known as 
hotspots. The overprcssured hotspot plasma then flows 
backwards, inflating a cocoon of radio-emitting plasma. 
Kaiser & Alexander (1997) argued that the pressure from 
the cocoon may be responsible for keeping AGN jets col- 
limated over long distances. These authors showed that, 
in this case, the radio source will necessarily expand self- 
similarly, its size growing as 



D = 2cii?o 
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(1) 



where ci is a constant of order unity (see Appendix), 
Qjet is the jet kinetic power; and the source is assumed 
to expand into an atmosphere with density profile p{r) = 

The observed radio luminosity is produced by syn- 
chrotron emitting electrons and/or positrons in the co- 
coon. Thus, the luminosity depends on the cocoon mag- 
netic field energy density ub, the distribution 71(7) of the 
Lorentz factors of the emitting electrons, and cocoon vol- 
ume V. In the absence of loss processes, it can be shown 
(e.g. Longair 1994) that, for an initial electron energy 

distribution ^(7^) = ke^~^ injected at the hotspot be- 
tween a minimum and maximum Lorentz factor 71 and 
72 , under the assumption of equipartition between par- 
ticle and magnetic field energy density, the cocoon radio 
luminosity at frequency v is 



5 + p I 



■{i+^r 



(2) 



where Ai is a numerical constant given by Equation I A2I 
in the Appendix. For our fiducial values of p = 2.5, 
71 = 10^ and 72 = 10^ we have Ai = 3.5 x 10"**. 



For self-similar expansion, the cocoon volume V — 
ttR^D^ where Rt is the axial ratio of the source. Again, 
with the assumption of equipartition between particle 
and magnetic field energies as the radio source expands, 
the magnetic field energy density is proportional to co- 
coon pressure, ub = Pc (r+i)(r -1) '^here r = ^^ and 
F/ is the a d iabatic index of the lobes (jKaiser fc Bestl 
I20071I2008D . IKaiser fc Alexandeil (|1997D give an expres- 
sion for the evolution of c ocoon pressure Pc wit h time 
(their Equation 20, see also IKaiser fc Bestl (j2007[) Equa- 
tion A4), which yields 



Ub 
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(r + l)(F,-l) 



where /p is a con stant given in Equation A5 of 
IKaiser fc Best! (|2007l ) (see also the Appendix). The radio 
luminosity is then 



L^^AittR^ 



rfp 



ir + l){Ti-l) 



'Pr)2 



X (1 + z)— [poR'oQtet] D 



,(l-p)/2 



3-(l+l)(^) 



(4) 



Note that the synchrotron spectral index is just a — 
(1 — p)/2. A reasonable hotspot spectral index between 
a = —0.5 and —0.75 gives p = 2.0 — 2.5. Because of 
radiative cooling discussed below, the observed cocoon 
spectral index is typically steeper. Taking p = 2 and 
assuming a reasonable value of /3 = 1.5 (e.g. Vikhlinin 
et al. 2006)0 yields 
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For this set of assumptions, at a fixed radio luminos- 
ity, the relationship between jet power and size is then 
Qjet oc D5/28_ Ag Willott et al. (1999) point out, for 
realistic values of D (typically 1 — 1000 kpc) this is a 

very weak dependence. Hence, one expects Qjot oc LJ^^- 
for powerful radio galaxies. We note that, as Equation [4| 
clearly shows, the Qjot — ^radio relation depends on the 
source spectral index. 

2.2. Size dependence of luminosity 

The above analysis is only correct if losses in the 
synchrotron-emitting population are negligible. Three 
loss mechanisms conspire to modify the electron energy 
spectrum: adiabatic expansion of the volume element in 
which the electrons were last accelerated, synchrotron 
radiation losses, and upscattering of Cosmic Microwave 
Background photons via the inverse Compton process. 
As a result of these various loss mechanisms, the Lorentz 
factor evolves as (c/ Kaiser et al. 1997 Equation 4) 
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where a2 = f (|±f 
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the energy density of the CMB photons, and the other 

^ A less conservative range of 1 < /3 < 2.5, 2 < p < 2.7 yields a 
range of [-1.2, 0.1] for the exponent of D in Equation[4l 
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Figure 1. Comparison of jot power estimates for the hotspot sam- 
ple. Hotspot jet powers are determined for g = 2 (see Paper 1), and 
are independent of radio source environment and orientation; while 
dynamical models assume the same environment, shape and orien- 
tation for all sources. These assumptions contribute to the scatter. 
Model parameters are as follows: cocoon axial ratio i?T = 2.3, the 
median value in our sample; minimum Lorentz factor 7niin = 1000 
IGodfrev et al.l 120091 ); and particle injection index p = 2.5. The 
density profile is as in Willott et al. (1999), with /3 = 1.5 and par- 
ticle density at 100 kpc nioo = 3000 m~'^. Shaded region shows 
the best fit between these estimates. The solid line indicates the 
expected 1;1 correspondence. 

symbols have their usual meanings. 

To a good approximation, the observed radio luminos- 
ity comes from electrons emitting synchrotron radiation 
at v = j'^h'L, where vl is the Larmor frequency. The 
evolution of 7(<) for a population of electrons with ini- 
tial Lorentz factor 7^ at acceleration time ti is then given 
by Equation [6l As shown by Kaiser et al. (1997), in- 
tegrating over all injection times then yields the radio 
luminosity as a function of source age. An important 
point is that the radio luminosity will evolve with age, 
as electrons radiate away their energy. This is not taken 
into account in the simple analysis of Section 12.11 

3. RESULTS 
3.1. Comparing jet power estimates 

Figure [T] shows a comparison of hotspot and dynamical 
jet powers for the sample described in Paper 1. These 
objects were selected from a complete z < 1 sample of 
FR II radio sources (jMullin et al.ll2008[) on the basis of 
having at least one (and often two) well-resolved terminal 
hotspots, allowing their jet powers to be determined. For 
details of the sample selection procedure and jet power 
calculations, we refer the reader to Paper 1. 

Dynamical jet powers were determined from source size 
and luminosity using a grid of dynamical models. These 
are generated as described in Shabala et al. (2008). 
Briefly, for each jet power, a source size and luminosity 
are calculated for each point in time given an environ- 
ment into which the radio source is expanding. Once a 
model grid is generated, the observed sizes and luminosi- 
ties are used to infer dynamical jet powers. Our modeling 
assumes all sources expand into the same environment, 
and have the same shape (i.e. axial ratio Rt relating 
cocoon length to its width). Formal uncertainties in dy- 
namical jet power estimates are largely due to uncer- 
tainty in the source orientation to the observer's line of 



sight. Error bars on hotspot jet power estimates corre- 
spond to the difference in derived jet powers in the two 
hotspots (see Paper 1 for details). While axial ratio val- 
ues arc available for individual sources, we chose not to 
use these in order to assess the sensitivity of dynamical 
jet power estimates to this parameter. Together with the 
"same environment" assumption, this yields some scat- 
ter in the Qjet,dyn — Qjet.hs relation. However, there is a 
strong correlation consistent with a linear relationship. 

3.2. Size-luminosity tracks 

In the following results we further restrict our sam- 
ple by requiring that jet powers from each hotspot pair 
agree within a factor of two. This ensures that the mea- 
sured hotspot jet powers are closely related to total ra- 
dio source size and luminosity, rather than being de- 
pendent on small-scale ISM inhomogeneities in or near 
the hotspot region. Our sample comprises of 44 radio 
sources. 

Figure [2] shows the size-radio luminosity tracks for this 
sample. We use hotspot jet powers to split our sample 
into high jet power (left panel), medium (middle) and 
low jet power (right) bins. Solid lines show the limit- 
ing expected size-luminosity tracks. For each panel, the 
lower track is obtained by taking the lowest jet power 
source in that bin, and evolving it at the highest redshift 
in the sample] (z = 1). Conversely, each upper track is 
obtained by evolving the highest jet power source at the 
lowest redshift (z — 0). 

Two things arc apparent from Figure [2] First, at a 
fixed jet power there is a clear dependence of radio lumi- 
nosity on source size. The shaded regions show radio lu- 
minosities expected from the size-independent relation of 
Willott et al. (1999; their Equation 12). Second, dynam- 
ical models do a very good job of reproducing this evolu- 
tion of radio luminosity with size. It is worth emphasising 
that the only parameter allowed to change between the 
three panels in Figure[2]is jet power, which is determined 
from a completely orthogonal set of (hotspot) observa- 
tions. All other parameters, such as those describing the 
density profile of the atmosphere into which the radio 
source is expanding, radio source shape, and initial par- 
ticle energy distributions, are held fixed. The decrease in 
radio luminosity with size is easily understood in terms 
of loss processes described by Equation \6\ 



3.3 



relation 



Mechanical power - radio luminosity - 

We can combine the three panels of Figure [2] to eluci- 
date the dependence of jet power on source size. Com- 
bining expressions for source size (Equation [1]) and radio 
luminosity (Equation |4|) yields 



-P 



logLi-adio-^rlogQjct = M- 



b-P 3(1 +p) 



P 



logD 



"(7) 

where A2 is a normalisation constant given in the Ap- 
pendix. The left-hand side of this equation represents 

■^ The redshift dependence comes from the conversion between 
the observed and emitted frequency, i^oba = i^em/(l + z). As radio 
lobes are steep-spectrum sources, the radio luminosity at a given 
observing frequency decreases with redshift. There is also a redshift 
dependence due to the changing CMB energy density (uc oc (1 H- 
2)*) in Equation |6] 
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Figure 2. Size-luminosity distribution of FR-II radio sources with 
well-defined hotspots. Jet powers are in the range 6 X 10'^*-3 X 10'^® 
(left), 6 X lO^'^-e X 10^* (middle), and 6 X 10^^-6 X 10^'' W (right 
panel). In each case, model tracks (solid lines) show the predicted 
evolution of highest and lowest jet power sources. For all plots, 
the parameters are as in Figure [l] Only the jet power is varied 
between panels. Predictions of the Willott ct al. (1999) relation 
are shown with a shaded region. 

a radio luminosity scaled by jet power in a ecordance 
with the dynamical model of lKaiser fc Alexande r (1997,). 
For our adopted parameters, Equation [7] reduces to 
0.533 log Lradio — | log Qjct = A2 ~ 0.233 logD. Because 
our sample spans a large redshift range {z ~ — 1), the 
radio luminosity must be adjusted accordingly. This is 
difficult to do analytically, but can be readily done with 
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Figure 3. Size-luminosity tracks of Figure [2] adjusted for jet 
power and redshift (see text). Only sources with Qj^t > 10** ergs 
s~^ are considered in this analysis. The IWillott et al] l|1999l ) re- 
lation predicts a horizontal line, ruled out by our Monte Carlo 
simulations at the 1 percent level. Medians are shown by open 
diamonds, and the best-fit relation of Equation [8] by the solid line. 

dynamical models. For our parameters, this yields a cor- 
rection of the forn0 Lradiol^) = ^radio,z=o(l + 2)^-^^- 

Figure [3] shows the dependence of the radio luminosity 
- jet power combination of Equation [7] on source size. 
We note that if loss processes associated with the radi- 
ating particles w ere not impo rtant, as is the case in the 
relation given by IWillott et a l. (1999,), we would expect 
all points to lie on a horizontal line. In other words, 
there would be no size dependence of the ordinate. We 
performed 10,000 Monte-Carlo bootstrap resampling re- 
alizations of our data, adopting an average uncertainty of 
(TiogQ^j — 0-2 dex in individual jet power measurements. 
This value comes from the observed distribution of north- 
to-south hotspot jet power ratios, plotted as error bars in 
Figure [3l Our Monte-Carlo simulations yielded a slope 
of -0.38 ±0.11 and intercept of 1.43 ± 0.28. The derived 
slope is significantly different from zero at the 1 percent 
level, ruling out the null hypothesis of no size depen- 
dence. As we explain in the following section, this result 
can be readily understood by considering the energy loss 
processes of the radiating particles. 

We use our data to provide a new estimator of jet power 
from radio luminosity and source size. 
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where we have assumed a hotspot jet power normali- 
sation constant g = 2, consistent with our findings in 
Paper 1. 

As we discuss below, the size term in Equation [5] plays 
an important role in the conversion of radio luminosity 
to jet power. In particular, flux and volume-limited sam- 
ples can yield very different scaling relations, due to the 
different correlations between source size and radio lumi- 

* This is the correction required to match the z = and z = 1 
size-luminosity tracks for the range of jet powers in our sample. 
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nosity inherent in each sample. 

4. DISCUSSION 

4.1. Spectral ageing 

The size dependence of the observed Qjet — iradio re- 
lation can be understood by considering the predicted 
relation in the lossless case (Equation [5|). Re- writing 
this in terms of the spectral index a = (1 — p)/2 gives 

iradio oc QJ^^""^/^D"('^)+(^-^). This is similar to the 
expression of O'Sullivan et al. (2011), but also includes 
a size term. Rearranging, 



Cjct - 
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where A3 is a function that absorbs the other model pa- 
rameters and numerical constants and is given in Equa- 
tion IA5I of the Appendix. As the radio source ages, 
loss processes discussed in Section 12.21 will cause the ob- 
served lobe spectral index aiobc to steepen to a value less 
than the injected spectral index ainj (note that "nor- 
mal" spectral indices are negative in our notation). This 
is a well-known observational result: for example, Blun- 
dell et al. (1999) find strong evidence for large sources 
having steeper spectral indices (see their Fig. 8), with 
sources larger than 100 kpc often having aiobc < —1.0, 
compared to aiobc ~ —0.5 for sub-100 kpc sources. This 
will affect the Qjot — ^radio relation in two ways. First, 
the steeper spectral index will affect the luminosity term, 
changing the slope from 3/(3 — a) = 0.87 (for a = —0.5) 
to 3/(3 — a) = 0.75 (for a = —1). Crucially, there will 
also be an additional size dependence. As argued by 
Willott et al. (1999), and in Section [^TTl this dependence 
is very weak for sources that have not suffered significant 
losses (aiobc ~ ainj ~ —0.5), with the exponent of D in 
Equation [9] ~ 0.2. Steepening of the spectral index to 
a = —1 makes the exponent of D equal to 0.5, consistent 
with the dependence found in Equation |S1 Importantly, 
for such sources an order of magnitude change in size 
(e.g. from 30 to 300 kpc) changes the normalisation of 
the Qjot — iradio relation by a factor of 3, with larger 
sources having higher jet powers at a given lobe radio 
luminosity. 

4.2. Effects of sample selection on the mechanical power 
- radio luminosity relation 

We illustrate the importance of the size dependence of 
jet power by considering another sample of radio AGN, 
a low-redshift (0.03 < z < 0.1) volume-limited sample of 
Shabalaetal. (2008). Figureljshows the relationship be- 
tween the dynamical jet power and 1.4 GHz cocoon radio 
luminosity. The sample is split in source size. The largest 
sources clearly lie above the best-fit relation, shown by 
the solid line. This is as expected from source dynamics: 
the oldest (and therefore largest) sources will suffer sig- 
nificant losses, and have a lower radio luminosity for a 
given jet power, than their more compact counterparts. 

We also plot the Qj^t — iradio relation for the hotspot 
sample of Paper 1 (blue dashed line). The two samples 
(hotspot and volume limited) clearly show very different 
scalings between radio luminosity and jet power, with 
slopes of 0.67 ± 0.05 and 1.02 ± 0.02 respectively. As we 
show in the top panel of Figure [5l the hotspot relation- 
ship derived in Paper 1 provides an excellent fit to the 
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Figure 4. Dynamical jot power - radio luminosity relation for the 
volume-limited sample of Shabala ct al. (2008). Despite significant 
scatter induced by environment and orientation effects, at a fixed 
radio luminosity the largest sources have systematically higher jet 
powers than their more compact counterparts. This systematic 
shift is primarily due to inverse Compton ageing of the electron 
population. The best fit Qj^t — ^1.4 relation for all points in this 
sample is shown by the solid black line. The dashed blue line is the 
relation derived in Paper 1 for the hotspot sample, consisting pre- 
dominantly of large (> 100 kpc) sources. The slope of the best-fit 
relation for the volume-limited sample is 1.02±0.02, compared with 
0.67 ± 0.05 for the hotspot sample. This difference comes about 
because of the different size distributions of the sources making up 
the two samples (Figure [Sjl. 

largest (>200 kpc) sources in the volume-limited sam- 
ple. This is because the hotspot sample consists mostly 
of sources of this size (see bottom panel of Figure [5]) . 
On the other hand, jet powers for the smallest, faintest 
sources in the volume-limited sample are overpredicted 
by 0.5 dex. 

The difference in the slopes of the Qjct — ^radio relations 
can be understood by considering the effects of the size- 
luminosity correlation in each sample on Equation [51 In 
the hotspot sample, the brightest sources are also the 
most compact (Figure [2|) . On the other hand, in the 
volume-limited sample it is the largest sources that are 
brightest (see Figure 12 in Shabala et al. 2008). Thus, for 
the hotspot sample the exponent of L151 in Equation [S] 
is < 0.8 while for the volume-limited sample it is > 0.8. 
The exact slope depends on the relative numbers of large 
and small sources, i.e. the size distribution shown in 
Figurc[5l Thus, any derived Qjot— ^radio relation depends 
sensitively on sample selection effects. 

We note that in this analysis we excluded all sources 
smaller than 10 kpc from our volume-limited sample. 
This is because in dense galaxy cores, source luminos- 
ity can actually increase as the source grows (Alexander 
2000). Hence for the smallest (< 10 kpc) sources, low ra- 
dio luminosity does not necessarily mean low jet power, 
and our dynamical model will systematically underesti- 
mate the jet power. It is worth noting that the number 
counts of these small sources (Figure [5]) far exceed those 
of larger AGN, arguing strongly in favour of a population 
that never grows to large sizes. These sources arc most 
likely disrupted due to a combination of low jet power 
and high environment density (Alexander 2000), and will 
exhibit a higher radio luminosity at a given jet power 
(sec Equation [9|), due to the so-called "environmental 
boosting" (Barthcl & Arnaud 1996). Thus, even with a 
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Figure 5. Distribution of dynamical jet powers for the volume- 
limited sample of Shabala et al. 2008 (left panel). The best-fit 
Qjet — iradio relation from the hotspot sample has been subtracted. 
The largest sources have dynamical powers broadly consistent with 
hotspot values (i.e. the black curve peaks near zero). On the other 
hand, for a given radio luminosity, smaller sources (red and purple 
triangles, and blue circles) arc seen to lie at systematically lower 
jet powers. The best fit is different for the hotspot and volume- 
limited samples. This is because the slope of the jet power - radio 
luminosity relation depends on the distribution of source sizes in 
the sample (right panel). 

proper treatment of their evolution, compact sources will 
have systematically lower Qjot/-^radio ratio, resulting in 
a further steepening of the relation in comparison with 
samples consisting of only large, powerful radio sources. 

4.3. AGN kinetic power budget 

We conclude with a brief discussion of the implica- 
tions of our findings on the AGN energy budget, and 
feedback from AGN on galaxies. Adopting a simple scal- 
ing relation between radio luminosity and jet power that 
docs not include a measure of source size or age results 
in significant errors in jet power estimates. For exam- 
ple, Equation [8] shows that calibrating the jet power - 
radio luminosity relation using the largest (> 500 kpc) 
sources results in overpredicting the jet power for the 
smallest (10 kpc) sources by an order of magnitude. 



Conversely, calibration of the relation using the small- 
est sources would underpredict jet power for the largest, 
oldest AGN by a similar amount. 

This is important, because understanding the detailed 
interplay between AGN jets and their surroundings is 
crucial to quantifying AGN feedback. Incorrect esti- 
mates of AGN jet power will result in mis-modelling of 
AGN feedback processes. For example, the most power- 
ful radio sources can affect multiple galaxies in a single 
AGN outburst (Rawlings & Jarvis 2004; Shabala et al. 
2011). Accurate jet power measurements, and interac- 
tion of jets with their environment, are crucial to estab- 
lishing the population of sources which can grow to large 
enough sizes to affect this feedback. On the other hand, 
a number of authors (e.g. Best et al. 2006; Shabala & 
Alexander 2009; Pope et al. 2012; McNamara & Nulscn 
2012) have argued that frequent low-level feedback from 
low-luminosity AGN is responsible for preventing excess 
cooling and star formation in the low-redshift Universe. 
Accurate jet power measurements in these objects are re- 
quired to test whether the heating and cooling rates are 
indeed balanced. Overall, our results indicate that accu- 
rate estimates of the integrated kinetic power output of 
AGN can only be obtained if a measure of radio source 
ages, such as size or spectral index, is used in addition 
to their luminosities. 

5. CONCLUSIONS 

We have used new kinetic jet power measurements 
of Godfrey & Shabala (2012) to examine the rela- 
tionship between integrated radio luminosity and jet 
power in steep-spectrum extragalactic radio sources. We 
find that the use of radio luminosity as a proxy for 
jet power under/over-predicts the jet powers for the 
largest/smallest sources. We provide a new fitting func- 
tion for jet power as a function of both source size and 
luminosity. 
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Because of this source size effect, any derived radio 
luminosity - jet power scaling relation is sensitive to 
the distribution of source sizes making up the sample. 
We have illustrated this point by comparing the sample 
of 3C FR-II radio galaxies used by Godfrey & Shabala 
(2012) with the complete volume-limited sample of Sha- 
bala et al. (2008). In the hotspot sample, the brightest 
sources are the most compact; the opposite is true for 
the volume-limited sample, where the largest FR-II radio 
sources are most luminous. The different correlations of 
source size with luminosity in the two samples yield very 
different slopes for the Qjct — -^radio relation, 0.67 ± 0.05 
and 1.02±0.02 respectively. Thus, caution must be taken 
when using radio source luminosity to infer AGN jet ki- 
netic power. 



Radio luminosity - mechanical power correlation 



APPENDIX 

NUMERICAL CONSTANTS 

We provide expressions for t he normalisat i on con stan ts discussed thr o ughou t the paper. These constants were given 
in iKaiser fc Alexander! (I1997D , iKaiser et al.l ()1997D and iKaiser fc BestI ()2007D , but for completeness we recount them 
here. We define F; and Fx to be the adiabatic index of the lobe and external medium respectively, p the energy index 
of the electron energy distribution n{"/) ~ kej~'^, Rt the axial ratio of the source, and the source is assumed to expand 

into an atmosphere with density profile p{r) = Po ( 7^ ) 

For models used in the paper, we assume values of Rt = 2.3, p = 2.5, /3 = 1.5, Fx = 5/3, Fj = 4/3, 71 = 1, 72 = 10^, 
p = 151 MHz. 
The constant ci in Equation [T] is 



ci 



i?rp 



(Fx + l)(Fz-l) (5-/3)3 



IStt, 



Ti + {Ti - 1) 



-4-/3 



1/(5-0) 



(Al) 



and ci = 1.08 for our assumed parameters. 
The constant Ai in Equations [2] and H] is 
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IG-K^Te 



(p+i) 



(2/io)" 
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f{pnui2)' 



(A2) 



where re is the classical electron radius, q and nie are the electron charge and mass respectively, fio is the permeability 
of free space, /(p, 71,72) = / ^ l^^^d'f, and the function 



C2{P)- 



3P/2 



Tfn (^ 



)rfn(f 



12 



rf„(f 



12^ 



p+13 p + 2 
2 2 TT 2 



= 2.26 X 10"^ for p = 2.5 

where Tfniz) = J^ t'~^e~*dt is the Gamma function (see e. g. iWorrall Il2009| ). For our parameters, Ai = 1.0 x 10" 
The constant fp in Equation [3] and SI which we quote directly from Equation A4 of IKaiser fc Best! (|2007l ). is 
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18^2(5-/5)73 



(Fx + l)(5-/?)2i?^ 



The constant A2 in Equation [7] is 

4 



A2 = 



log 



AittR^ 



rfp 



(r + l)(Fz-l) 



PaRo 



(A3) 



(A4) 



where r ^ {p + l)/4. 
For our choice of parameters, A2 



-5.9. We note that this expression applies only to the lossless model of radio 



source evolution, and this value will decrease for realistic evolutionary models. 
The constant A3 in Equation [9] is 



A3 = {AittR^) 



rfp 



(r + l)(F,-l) 



PoRq 



(A5) 



where A3 — 8.8 x lO'^ for our choice of parameters. The value of A3 will increase when losses are taken into account. 
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